Testing DM halos using rotation curves and lensing: 
A warning on the determination of the halo mass. 
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There are two observations of galaxies that can offer some insight into the nature of the dark matter (DM), 
namely the rotation curves and the gravitational lensing. While the first one can be studied using the Newtonian 
limit, the second one is completely relativistic. Each one separately can not determine the nature of DM, but both 
together give us key information about this open problem. In this work we use a static and spherically symmetric 
metric to model the DM halo in a galaxy or in a galaxy cluster. The metric contains two free functions, one 
associated with the distribution of mass and the other one with the gravitational potential. We use galactic, 
typical rotation curves to univocally determine the kinematics of the halos. We compute separately the mass 
functions for a perfect fluid and a scalar field, and demonstrate that both models can be fitted to the observations, 
though with different masses. We then employ lensing to discriminate between these models. This procedure 
represents a test of models using two measurements: rotation curves and lensing. 



o 

Oh 
I 

o 



> 

in 

00 

o 
o 



There are at least four main groups of cosmological ob- 
servations which are consistently described by the ACDM 
model, namely: the acoustic peaks in the cosmic microwave 
background radiation, Supernovae type la data, rotation 
curves of spirals and dynamics of galactic clusters, and their 
gravitational lensing (li|-[3|]. The last two of these observa- 
tions are related most significantly with the presence of dark 
matter (DM), where a consistent model considers a DM halo 
surrounding the galaxies and galactic clusters. 

Using Newtonian gravity we can infer the mass of the DM 
halos from observations of velocity profiles in galaxies and 
galaxies clusters. Alternatively, one can use lensing images 
or distortions of background galaxies due to the space-time 
curvature iHHtl- This is notorious, particularly in the case of 
gravitational lensing, which is a general relativistic effect per 
excellence, as it strictly can not be dealt with using Newtonian 
physics. 

There has been many attempts to model DM halos in the 
general relativistic framework, some of them has shown that 
general relativistic effects can be important ll6l- [l3ll . but in gen- 
eral they make use of approximations or limits that in some 
cases represents assumptions about DM properties. The scalar 
field DM model is an alternative proposed in the past 1 14] to 
fit the observed amount of substructure Iil5i] . the critical mass 
of galaxies 1 16], the rotation curves of g alaxies iUtI] . the cen- 
tral density profile of LSB galaxies [18], the evolution of the 
cosmological densities [19.1 , among other topics. 

In the present work we focus on two models for the DM 
halo (a perfect fluid and a scalar field) and explicitly show 
what exactly is determined by the observations and what 
comes as extra assumptions. This is particularly important 
as the unknown nature of the DM is one of the most rele- 
vant questions that one would like to solve. It is clear that to 
make assumptions specifically on the nature of DM, in order 
to obtain information on the nature of the DM, is skewing the 
problem. 

In our study we consider a static and spherically symmetric 
space-time in General Relativity, described by the line ele- 



ment: 
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where dQ. - d0^ -H sin^ Odip^. The gravitational potential 0(r) 
and the mass function m(r) are functions of the radial coordi- 
nate only. In fact, due to the symmetries of this space-time, all 
physical quantities depend only on r. The Einstein's equations 
are the known set [|20ll : 
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where prime ' = didr. The potentials of such space-time are 
determined by the DM halo. In order to see the difference 
between two hypothesis on the nature of DM, we will consider 
two types of composition for halos: a perfect fluid and a scalar 
field. 

The above equations are complemented by the conservation 
equation of the matter-energy generating the curvature of the 
space-time. But given the different nature of the fluids consid- 
ered here, this equation is treated separately for each fluid. 

In the case of the perfect fluid, the stress-energy tensor is 
given by r^y = {p -\- plc^)u^Uy -n pg^v, where the density is 
p - {\ -\- e)po, where po is the rest mass energy density and 
€ the internal energy per unit mass, is the co-movil four 
velocity, normalized as m^m'^ - -c^, and p is the pressure. The 
conservation equation, T'^y;^ = 0, implies the field equation: 
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which can be rewritten as 
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Now, for the scalar field the stress energy tensor is given by 
1 
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where - d(p/dx"; and V ((p) is the scalar potential. The 
components of the stress-energy tensor are 
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From the conservation equation for the scalar field, T'^y-f, 
0, one obtains a field equation, the Klein-Gordon equation. 
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that can be written as 
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which is remarkable similar to the field equation for the per- 
fect fluid, Eq.©. Given this similarity, it is convenient for our 
mathematical description to consider the single field equation 
for both types of matter 
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in which a - for the perfect fluid, and a = 1 for the scalar 
field. Notice that if one considers a sort of perfect fluid given 
by T^y = diag(-p, p, pt, pi), pi (some times called "tangen- 
tial" pressure) denotes a term representing the ignorance we 
have on the features of the fluid. This presure pi is related 
to the other fluid variables as pi - (1 - a) p - apc^, (see 
Eq. ((Toll), where "a" takes, in principle, any value. There are 
works which have discussed this field equation considering a 
as a free parameter ifTTUTsll . For the purpose of the present 
work we will consider only the two extremal cases, a = 
and a = 1, but the discussion can be directly applied for these 
cases as well. 

In this way, the system of equations which must be solved 
are, the Einstein's equations, Eqs. (|2|3t , and the field equation, 
Eq. ([Tol l. In either case, there are four unknown functions, 
m, <1), p and p, for the case where the curvature of the space- 
time is due to the perfect fluid, or m, O, </> and V((f)) when the 
curvature is caused by the scalar field. Thus, we have three 
equations for four unknown functions. In either case, we need 
only one extra data. It is important to underline this fact. Once 
the extra data is given, there is no more room left for any other 
assumption, the rest of the functions are determined by the 
system of equations. If, for instance, we give an equation of 
state for the perfect fluid, p = p(p) or, in the case of the scalar 
field, an explicit form for the potential, V((p), there is no free- 
dom left to choose the form of the rest of the functions, they 
will be determined by the system of equations. 

Following the line of work presented in |8], we use obser- 
vational results to close the system of equations. In the case 
of galactic halos, two main observations can serve to obtain 
the desired information: measurements of rotation curves in 



spirals and light deflection by lensing. In this work we choose 
the former to complement the above field equations and use 
the latter to discriminate between dififerent halo type models. 

Rotation curves. The motion of test particles in such space- 
time is determined by the geodesic equations and, for test par- 
ticles in circular motion, there is a relationship between the 
gravitational potential, O, and the tangential velocity of those 
particles, Vc'. 
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where we have definedyS^ = This tangential velocity is the 
one measured by observations of rotation curves in galaxies. 
Thus, Vc is an observable function, and by means of Eq. ( fTTT ). 
the gravitational potential can be determined. Thus, given this 
observable, there is no room left for an equation of state for 
the perfect fluid or for a given scalar field potential. 

Moreover, as long as the magnitude of the observed ve- 
locities are small with respect to the speed of light, this jus- 
tifies the validity of one of the weak field approximations 
O/c^ << 1 that one usually assumes by taking the weak field 
limit. Here we want to emphasize that the approximations 
2Gm I r « 1 and especially p « p are, in general, extra 
hypothesis which strongly depend upon the nature of the DM 
type. It is clear that if all these conditions are satisfied, then 
the above system of equations, Eqs. (|2][3lllll, reduces to the 
hydrodynamic set of equations for the case of the perfect fluid 
model. But, for example in the case of the scalar field, there 
is no Newtonian limit, and one has to be careful with these 
approximations. 

After substituting Eq. ( fTTT i into the gravity equations, Eqs. 
(HO, we obtain an equation (with no approximations) for the 
mass function as the only free function 
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The functions P(r) and Q(r) depend on the type of fluid we 
are dealing with (a) and on the rotation curves profile. 

The mass function can be expressed in terms of the gravita- 
tional potential, O, through the integral 



J'P(r')dr' 
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where the value of the integration constant, C, is set by the 
appropriate boundary conditions. 

For the case of the perfect fluid, the density and pressure 
are directly computed from Eqs. (|2]i and (O, respectively. 
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For the scalar field, using the expressions Eqs. (Q, we ob- 
tain that 
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Once the function j6(r) is given, the scalar field and its po- 
tential are straightforwardly determined in terms of the radial 
coordinate. In order to obtain the form y(0), one needs to 
invert the solution for the scalar field (r - r{<p)), and to substi- 
tute it into Eq. (fTSI l. As shown below, this procedure works at 
least for simple /3 functions. 

In this way, we have shown that the mass function m, as- 
sociated to a galactic halo by means of the rotation velocity 
strongly depends on the DM model that is being considered. 
The single observation of the rotation curve is not sufficient 
to determine the nature of DM and hence the mass associated 
with the halo. Moreover, we have shown that the relationship 
between the pressure and the density, or between the scalar 
field and the scalar potential is fixed, up to integration con- 
stants, once the rotation velocity is employed. 

Lensing. The other observation concerns the gravitational 
lensing, that for the line element given by Eq. ([TJ, the deflec- 
tion of the light ray, A ip, at the radius of maximal approach, 
r,„, is given by ll2ll] . 
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Since the gravitational potentials and the fluid variables are 
already determined by the rotation curves of spirals, deflection 
angle measurements can serve to discriminate between mod- 
els. Here we deal with two examples, perfect fluid and scalar 
field DM models. Yet, observations of spirals that lenses light 
are not very common, however the first examples of them has 
recently appeared f23\. We remark that from the expression 
of the deflection angle, it is a large step to infer the mass of the 
DM halo based solely on the observation of the deflection an- 
gle. A supposition has to be made on the relation between the 
gravitational potential, O, and the mass function, m ll2lll . Such 
supposition, as we have shown, not only strongly depends on 
the type of matter, but also on the specific characteristics of 
the type of matter considered. 

In the next section we present some examples of known ro- 
tation curve profiles to give a quantitative description to these 
conclusions. 

Examples. The idea in this section is to stress the conclu- 
sions that we are presenting by means of considering a typical 
observation of rotation velocities in spirals and to directly de- 
termine the gravitational mass in each case, when the DM is a 
perfect fluid (dust) and when it is a scalar field. 

In practice we can consider a velocity distribution, as a phe- 
nomenological model, for instance the velocity profile coming 



from N-body simulations given by NEW 11231 124|] or a Burkert 
profi le |[25i1 given by the phenomenological of rotation curves 
12a], to determine the mass of each type of fluid. We will 
show that the gravitational mass inferred by the same veloc- 
ity profile is strikingly different for the perfect fluid and scalar 
field cases. 

Constant velocity profile. We will consider the simplest 
case of constant rotation curves as our first example. Although 
there are some examples of galaxies that present a constant 
velocity profile, for a few disk length scales |27|], this is not 
a typical behavior, being our own Galaxy a good counter ex- 
ample 1 28] and, in fact, there is an important rotation curve 
phe nomenology described by the Universal Rotation Curve 
II26II29I1 . However, the constant velocity profile offers us the 
mathematical simplicity to obtain straightforward analytical 
results and to show the main point of our work. 

For the gravitational potential, from Eq. ( fTTT i. when the ve- 
locity function is a constant, /3q, we get 
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The mass function can be analytically obtained for any value 
of the parameter a as: 
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where C is the integration constant of Eq. ( fT4b . For the case of 
DM described by a perfect fluid, a = 0, we fix this constant to 
zero in order to avoid changes in the signature of the line ele- 
ment, Eq. ([T]i. Thus, the mass function, and the corresponding 
pressure and density in the case of the perfect fluid are given 
by 
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We can see in the limit of very small velocities, ySo << 1, we 
recover the Newtonian limit, and pressure is negligible with 
respect to the density as we mentioned above, however, it is 
not zero and, actually, we obtain a barotropic equation of state, 
p = Wop. 

On the other hand, considering the DM halo due to a scalar 
field, the mass function is obtained from Eq. (fT9l ). with a - I. 
In this case, the mass function has a very peculiar behavior 
The first term is small, proportional to y6o"^, but negative. The 
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second term, proportional to the constant C, goes as 
thus, by choosing a positive value for the constant C, one can 
have a positive mass function for a large region, but this func- 
tion will present a divergence at the origin. Of course this re- 
sult was expected, as the space-time metric is static. In order 
to avoid this problern, we had to take non-static space-times, 
like the oscillatons fl6\, but this is beyond the scope of this 
work. It can be shown, however, that this divergence is cov- 
ered by an apparent horizon. Some features of this case of 
scalar field with a non zero constant C in the mass function, 
have been discussed in [11]. For the purpose of this work, 
we only notice that the geometric functions and those of the 
scalar field, have a non intuitive behavior, but are consistent 
with the rotational curve. Explicitly, for the case of C = 0, the 
mass function is 



y8o 



(24) 



and, with the geometric functions determined, the scalar field 
and scalar potential are completely fixed, given by: 
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where the expression for the scalar field, 0(r), was inverted to 
obtain r(0), and then express the scalar potential in terms of 
0, as explained previously. The "effective mass" of the scalar 
field, Weff -J , depends inversely on the characteristic 

distance of the halo. This distance is of the order of kilo- 
parsecs, and A ~ 10""^, thus it will turn into a typical mass for 
scalar field in a galaxy, which corresponds a very light boson 
mass ~ 10"^^eV/c^. This result is in agreement with the one 
obtained in previous works, see for example 1 15]. 

Going back to our previous discussion, notice how remark- 
ably different are the mass expressions derived from each type 
of fluid, being both consistent with the observed rotation ve- 
locities. This is the simplest case in which we can show how 
the single observation of the rotation velocities in halos de- 
termines the features of the perfect fluid model or the scalar 
field. 

Although the mass associated to the scalar field results neg- 
ative and this can be taken as a no-go result for static scalar 
field halos llsoll . rotation curves of spirals are not exactly flat 
(see discussion in {tS, '29^) and, in addition, we have to be 
cautious with the supposition of a static metric which is very 
restrictive for the scalar field . Thus, the above-result should 
not be taken as definitive, at most, it should be taken as a re- 
mark that a static DM halo is not well described by a static 
scalar field. A negative mass, or positive gravitational poten- 
tial, is known since long time ago IsT] from the fits to rota- 
tion curves using modifications of newtonian gravity in which 
a scalar field induces a Yukawa-type force. At the end, de- 
manding a constant velocity profile all the way in the radial 



direction implies an effective repulsive force to be acting on 
test particles in the galaxy. 

In any case, it emphasizes our point in showing how 
strongly depends the determination of the mass of the DM 
halo on the type of matter considered to describe it. 

The deflection angle for the case of constant rotation veloc- 
ity, considering the perfect fluid and the scalar field with the 
constant C -Q, implies the following expressions: 
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where we have defined x - Since in any case the deflection 
angle is a constant, i.e. it does not depend of the maximal 
approach radio r,„, it can be evaluated for a typical value of 
the velocity. For comparison we take the value A) = 1/1200, 
that corresponds to a velocity of Vc - 250 km/s. Evaluating 
the deflection angle, we get 



apf = 0.899547, 
a,f = 0.449546, 



(31) 
(32) 



for both cases the deflection angle is given in arc second units. 
We can see that there is a difference of almost half arc sec- 
ond between them, and the simultaneous observation of the 
rotation velocity and the deflection of light produced by the 
galactic halo, can teach us about the true nature of the DM. 

We notice that the deflection angle for the scalar field with 
a non-zero value of the constant C in the mass function takes 
very large values, a fact which certainly allows us to discard 
this option as a model for the DM halo, independently of any 
interpretations of the mass function. 

Now we study an example that is less striking though. 

Navarw-Frenk-White (NFW) velocity profile. Indepen- 
dently from its origin, the NFW profile [23, 24] is considered 
per se as a viable fitting model to describe galactic kinematics. 
This profile has been subject to geometrical studies elsewhere 
lf32i[32l[33ll . In this example, we assume this profile as a valid 
phenomenological galactic profile for the galactic data. We 
obtain the usual expression for the mass derived within this 
description, and compare it with the same form of the rotation 
velocity, but considering that it is due to a DM halo composed 
of a scalar field. 

The velocity profile in the NFW model lEsi l24ll is given by 



r/rp 
1 -t- r/ro 



+ ln[l +r/ro] 



(33) 



where ctq - AkG po ro^ is a characteristic velocity of stars in 
the halo, given in terms of a characteristic density, and ro is 
a scale radius. Given this velocity profile, we have to solve 
Eq. ( fT4l i with a - Q for the perfect fluid and with a = 1 for the 
scalar field. In neither case there is an analytical solution, thus 
we have integrated the equations numerically. We do not want 
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to treat here specific galaxies but to emphasize the differences 
between the galaxy models. Therefore, we set ctq and ro to 
some typical values. In our plots we assume geometric units 
(G = c - 1), and therefore the characteristic velocity takes 
values, < ctq < 1, and the mass is less than the unity. For 
definiteness, we assume ctq - 0.001 and ro = 1. In figure 
([TJ we plot both halo masses (perfect fluid and scalar field). 
Disregarding the behavior near the origin, as long as we are 
considering the outside region, as mentioned above, we see 
that the mass associated to the halo in each case are differ- 
ent. We now consider lensing. By integrating Eq. (fTTT l for 
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FIG. 1 : Comparison of the masses using the rotation velocity profile 
from the NFW model with a perfect fluid (upper curve) and scalar 
field (lower curve). 

the given rotation velocity, Eq. ( |33] |, we obtain the following 
expression for the gravitational potential: 



-o-oro- 



(34) 



We substitute this expression, together with the corresponding 
numerical solution one for the mass in each case, in the equa- 
tion for the deflection angle, Eq. ( fTTb . and perform the integra- 
tion varying the value of the radius of maximal approach, r,„. 
The results are plotted in figure |2] As we see, the observation 
of the deflection angle can determine which type of matter is 
actually composing the DM halo. 



FIG. 2: Deflection angle generated by the gravitational lensing of a 
NFW rotation profile with a perfect fluid (upper curve) and scalar 
field (lower curve). 

In this example, the mass associated to the scalar field 
model is essentially positive and a well-behaved function, that 
is, it does not follow the no-go result mentioned in the previ- 
ous example. Though there is a small region where the mass 
becomes negative. This is due to the fact that we are demand- 
ing the velocity profile to grow in that region. In a real setting 
however the stellar disc adds to the velocity profile, thus we 
expect that its contribution avoids negative mass regions for 
the scalar field. 

With these examples it is clear how two different types of 
matter (perfect fluid and scalar field) can be consistent with 
the observation of rotation curves of DM halos, though they 
lead to different conclusions to the mass function inferred. 
The deflection of light can then be used to discriminate be- 
tween the two models. Even though the mass function for 
some model has not an intuitively expected behavior, it is nec- 
essary to use the observation in order to discard the model, 
being aware of the assumptions made during the derivation of 
such conclusions. 
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